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RESUMEN 

En este trabajo se presenta una comparación de tres métodos para deducir una propiedad de los 
potenciales de los elementos que permite exprear en distintas formas el potencial de Gibbs de los 
sistemas en equilibrio químico. Estos tres métodos son: (1) el procedimiento algebraico 
ampliamente usado en la modelización de las operaciones de destilación reactiva que lleva a la 
definición de las variables de composición transformadas, (2) la clásica formulación 
no-estequiométrica de los problemas de equilbrio químico, y (3) la transformación de Legendre 
completa del potencial de Gibbs con respecto a las variables de composición de los elementos. 
Estos tres métodos conducen al mismo resultado en la mayoría de los casos de interés práctico, es 
decir cuando los elementos se toman del conjunto de los componentes. El último, por su 
consistencia con otras formalizaciones de la Termodinámica, permite comprender mejor el 
significado de las variables de composición transformadas y de la función auxiliar del método de 
los multiplicadores de Lagrange. Las variables de composición de los elementos son las variables 
conjugadas de los potenciales de los elementos. La función auxiliar del método de los 
multiplicadores de Lagrange es la transformada de Legendre completa del potencial de Gibbs con 
respecto a las variables de composición de los elementos y su valor es cero en los estados de 
equilibrio. 
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1. INTRODUCCIÓN 

En las mezclas reactivas, los cambios en las variables de composición no son 
independientes. Por ello, para sistemas cerrados en equilibrio químico, la 
expresión diferencial de la ecuación fundamental adopta una forma más simple 
que en el caso de los sistemas abiertos no reactivos (Luetich 2002). Se ha 
demostrado que esta simplificación se debe al hecho de que esas variables pueden 
ser convertidas en otras que no cambian cuando el equilibrio es desplazado 
modificando las variables de estado. 

Un método algebraico para llevar a cabo esta transformación fue introducido por 
Barbosa y Doherty (1987) con el propósito de transferir a la destilación reactiva 
las técnicas de modelización de la destilación no-reactiva. 

Pérez Cisneros y otros (1997) trataron los mismos sistemas por medio de la 
clásica formulación no-estequiométrica (Smith y Missen 1991), mostrando que 
este método hace posible el uso de variables de composición que no avanzan 
incluso en el caso de reacciones tales como las de esterificación, es decir cuando 
no todos los elementos existen libres en la mezcla. 

Al mismo tiempo, Alberty (1997) presentó un método basado en el uso de 
transformaciones de Legendre para tratar los sistemas que contienen un 
componente cuyo potencial químico permanece constante. Su enfoque estaba 
basado en trabajos previos relacionados con la reducción del número de variables 
de composición en sistemas que contienen seudoisómeros (Alberty y Oppenheim 
1990). 

En este trabajo se demuestra que el uso de transformaciones completas de 
Legendre del potencial de Gibbs con respecto a las variables de composición de los 
elementos permite hacer una interpretación más profunda del sentido físico de 
las variables de composición transformadas (Luetich 1999) y de la función 
auxiliar del método de los multiplicadores de Lagrange (Luetich 2001). 

En el desarrollo se hace una comparación de las tres técnicas, aplicándolas a la 
solución de un problema en el cual conducen al mismo resultado. 

Este artículo es el segundo de un grupo de cuatro trabajos relacionados escritos 
para clarificar el concepto de reversibilidad microscópica. Con este objetivo, el 
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foco se pone aquí en la información contenida en las distintas formas de la 
ecuación fundamental para sistemas en equilibrio químico. 

2. ESCRITURA DE LA EXPRESIÓN DEL POTENCIAL DE GIBBS 

PLANTEO DEL PROBLEMA 

Dado el sistema reactivo estudiado por Brandani y DiGiacomo (1984) 



(2.1) 



i formaldehído + agua ^ metilenglicol 
v metilenglicol ^ agua + polioximetileno 

se trata de hallar la expresión más simple para el potencial de Gibbs. 

SOLUCIONES 

■ COMBINACIONES LINEALES 

Para simplificar la notación, se hará uso de componentes genéricos A, B, C, y D, 
en lugar del formaldehído, el agua, el metilenglicol y el polioximetileno, 
respectivamente. 

( A+B^C 

Uc^B + D (2 ' 2) 

En el equilibrio, el potencial de Gibbs del sistema (2. 2), 

G = (fi A He He y"D ) ( n A n B n c n B ) T , (2. 3) 

está sujeto a las restricciones 



1 1 


\ 


1 


-1 


-1 


2 


l o 


-J 



( yU A jU B Mc jUd ) , „ = ( ). (2. 4) 



(Esta propiedad justifica el uso del adjetivo "químico" para la magnitud ¡x.) 

Entonces, sólo dos potenciales químicos son independientes. 

(1 \ (-1 2 \ 

(yU A /íb)^ .J + Í^c f*v)[ q _ 1 J = (0 0) (2 ' 5) 

1 V-l 2 ^ _1 



(fie fio) = ~(Ha Mb)\ II _ | (2.6) 
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C y D son los llamados componentes de referencia (Barbosa y Doherty 1987). 

Por otro lado, asociando términos en la expresión matricial del potencial de 
Gibbs, 



m rn 

G = ( yU A Hb)(ji a n B ) + ( Me l¿T>)(n c n B ) . 



(2.7) 



El siguiente paso es el más importante debido a que, por medio de él, se eliminan 
dos variables: 



G = ( ¡x A A*b ) I 



n A + n c + 2n B 
n B + n c + n B 



(2.8) 



Esta relación permite definir las variables de composición transformadas (ñ) . 

ñ A \ a 



(Ha Hb) 



n B 



G 



(2.9) 



En el apéndice A se demuestra que este procedimiento es equivalente a eliminar 
los grados de avance de las dos reacciones. 

■ MULTIPLICADORES DE LAGRANGE 

El sistema puede ser visto como si estuviera constituido por los "elementos" A y B 
(aproximación física): 

f A + B ^ AB 
12AB^ B+A 2 B 



(2. 10) 



Por lo tanto, se trata de minimizar la función 

G = G(T,P,(n A n B n c n D )), 
sujeta a las condiciones 



(2. 11) 



n, 



n\ 





(n A \ 


1 1 2\ 


n B 


1 1 lj 


n c 




< n D j 



(2. 12) 



que definen las variables de composición n* A y n* B . Estas variables surgen de 
inspeccionar la cantidad de átomos que conforman cada molécula. 
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El problema es el clásico problema de minimización con variables ligadas que 
Lagrange resolvió haciendo uso de multiplicadores. 

Sea 



O = G- 



( A A A B ) 



10 12 
111 



ín K \ 
n B 
n c 

V n B ) 



(2. 13) 



la función auxiliar correspondiente [el lagrangiano, 0(n, A)]. 



Sus derivadas con respecto a los números de moles son 



( ^° 1 

dn A 








50 




(Ha) 




dn B 




Hb 




50 

dn c 

50 




He 

v Hv t 




V dn B ) 









- (A A A B ) 



10 12 
111 



(2. 14] 



El producto de matrices debe ser transpuesto porque 

d _ 

(n A \ 




n B 



(10 0) J 



(0 10 0) J 
(0 1 O) 1 
(0 l) 1 



(2. 15) 



Estas derivadas se hacen nulas — llevando G a su valor mínimo — cuando 

( A A ■ 

A B 

Aa + A b 

2 A a + A b 



(Ha) 




(1 0^ 
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Hb 
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í Aa 1- 


He 




1 1 


UbJ" 


vHn j 




12 lj 


V 



(2. 16) 



Por lo tanto, 



mín G = ( n A + n c + 2 n D n B + n c + nn) 



( n* A n* B ) | A I = mín G. 



Ap 



(2. 17) 

(2. 18) 
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En el apéndice B se muestra la relación entre los dos métodos precedentes. 



TRANSFORMADAS DE LEGENDRE 



Para este sistema reactivo, el potencial de Gibbs es: 



G = Oa n B n c n-o) 



V A<d ) 



(2. 19) 



A T y P constantes, la ecuación que determina el cambio G con las variables de 
composición es: 

n B 



(dG) TP = ( yU A jU B yU G yU D ) d 






(2. 20) 



Por lo tanto, el problema consiste en encontrar el número mínimo de veces que el 
potencial de Gibbs debe ser transformado — cuando las variables están ligadas 
por la estequiometría de las reacciones químicas — para obtener la transformada 
de Legendre completa. Pero este número es precisamente el número de 
elementos. En otras palabras, las variables que permanecen constantes durante 
una reacción son las variables de composición de los elementos. Por lo tanto, se 
debe hacer uso de la transformación (2. 12). 



En este caso, la transformada de Legendre completa de G es: 



G-{n\ ti*) 



ÍB 



G* 
i 



AB- 



(2.21) 



[Aquí se introdujo un nuevo símbolo (£) porque no se sabe todavía qué variables 
deben acompañar a las variables de composición transformadas. Las variables £ 
son los potenciales de los elementos.] 



Operando sobre las ecuaciones (2. 21) y (2. 12) 



d\G-{n\ ti* b ) 



6 



0; 



(2. 22) 



Actas Acad. Luv. 2002, 1, 23-36 



El potencial de Gibbs de mezclas reactivas 



29 



(fÍA-U /^B-ÍB y"C-ÍA-^B ^D-2^ A -Íb)C? 



n B 
n c 

\ n B ) 



= 0. 



(2. 23) 



Dado que esta condición se debe cumplir cuando no todas las variaciones son 
nulas, es posible obtener las expresiones para los potenciales de los elementos y 
compuestos: 



( ¡i A fi B ¡j, c fin ) = ( £ A & £a + £b 2 ¿f A + ¿b )• 



(2. 24) 



Por lo tanto, 



G = {n\ n\) 



Ha 
Hb 



(2. 25) 



Éste es el mismo resultado obtenido por el método de los multiplicadores de 
Lagrange. Nótese que este procedimiento — usando ecuaciones integrales — es el 
inverso del método de las combinaciones lineales (apéndice C). 

De aquí en más, para completar la deducción, será conveniente hacer uso de la 
siguiente notación: 

( Ha - U /¿b - & He ~ U - ft Hv - 2 ¿"a - & ) = ( H * A H*b H*c H*o ), ( 2 - 26) 
donde \x* es potencial transformado. 

Además, la transformación (2. 21) define un potencial termodinámico cuya 
expresión, en términos de las variables originales y transformadas, es: 



G*ab = ( ¡Je - Ha - I¿b Hn - 2 ¡j. a - Hb ) 



n B 

V n D ) 



(2. 27) 



G*^ = (0 n* c //* D ) 



n B 



(2. 28) 



Este potencial se hace cero en los estados de equilibrio. 



Entonces, la forma diferencial de esta función es: 



Actas Acad. Luv. 2002, 1, 23-36 



30 



J. J. Luetich 



(dCr ab^p 



, i df¿ A 1 



(2. 29) 



En la tabla 1 se muestran las fórmulas que relacionan los potenciales 
termodinámicos original y transformado y sus derivadas con las variables 
transformadas para los estados de equilibrio. 



Tabla 1 

¡ A+B^C 
Condiciones de equilibrio para el sistema < en términos de los potenciales de Gibbs 

^2L ^ d + D 



original y transformado 



G = mín G = n\ ¡x A + n* B ¡x B * 
(cIG)tp = fj, A dn* A + yu B dn * B = 



G* m = mín G*^ = n* c n c + ¡j,* d n D = 
(dG*^)^ = -n\ dn K - n* B d/u B = 



Ambos potenciales (G y G*^) tienen diferenciales nulos en los estados de 
equilibrio, pero sólo el potencial transformado es nulo. En el primer caso, el 
hecho de que el diferencial sea nulo se debe a la invariabilidad durante la 
reacción de las variables de composición transformadas. En el segundo caso, esto 
es así debido a que las variables transformadas están ligadas a los potenciales por 
una relación análoga a la ecuación de Gibbs y Duhem. 

■ Generalización 

La expresión general del potencial de Gibbs para sistemas en equilibrio químico 
es: 



G = G (m, n) = n T /i, 



(2. 30) 



donde m es el vector de las variables de estado y n es el vector de las variables de 
composición. 



m = (T, P) 



(2.31) 



Las restricciones estequiométricas reducen el número de variables de 
composición: 



Cn = n*. 



(2. 32) 



C es la matriz de conservación de los elementos. 
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Las variables intensivas que acompañan a las variables de composición 
transformadas en la expresión abreviada del potencial de Gibbs pueden ser 
relacionadas con los potenciales químicos. 

G = (n*) T £ (2. 33) 

V* = H-C T £=0 (2.34) 

Las formas diferencial e integral generales del potencial de Gibbs transformado 
son: 

G* = (jt *) T n; (2. 35) 

(dG*) m = -(n*) T o^. (2.36) 

De las ecuaciones (2.34), se obtienen las componentes de {. Cuando éstas se 
reemplazan en las ecuaciones (2.36), proporcionan la expresión de dG* en 
función de los potenciales químicos originales (//). La fórmula del potencial 
transformado resulta de las ecuaciones (2.34). Si ellas se reemplazan en la 
ecuación (2. 35), se obtiene la forma integral de G* . 

De otro modo, agregando a la ecuación (2. 23) la ecuación de Gibbs y Duhem, 

(li*) T dn + ndf¿ = 0; (2.37) 

[d (G - (n*) T ¿)] m = (dG*) m = - (n*) T d£. (2. 38) 

Por lo tanto, se obtiene la ecuación: 

(dG) m = fdn*, (2.39) 

cuya importancia se pone de manifiesto al resolver problemas de equilibrio 
químico en el marco de la aproximación física (Luetich 2002). 

Las variables de composición transformadas son las variables conjugadas de los 
potenciales de los elementos (Luetich 2001). 

La tabla 2 presenta un resumen de los tres métodos aplicados a la solución del 
problema planteado. 



Actas Acad. Luv. 2002, 1, 23-36 



32 



J. J. Luetich 



Tabla 2 

Un problema, tres enfoques 



ELPROBLEMA 


FORMAS ALTERNATIVAS 


Dado el sistema reactivo 


Combinaciones lineales 


, A+B^C 
120B + D 


Transformar el número de moles de C y D 


para expresarlos como combinaciones 


se trata de hallar 


de los de A y B. 


la expresión más simple 


G = n A fi A + n B ^i B + n c ¡i c + n D fi D 


para el potencial de Gibbs . 


í fJ-A + j"B - j"C = 
l -Hb + 2 fl C - yU D = 






Multiplicadores de Lagrange 




Usando las variables de composición 




de los elementos n* A y n* B , 




obtener la expresión del mínimo 




de la función 




G = G(T,P,(n A n B n c n D )). 




Transformadas de Legendre 




Encontrar las variables conjugadas 




de las variables de composición de los elementos 




para escribir la expresión diferencial de 




G = G(T,P,(n\ n* B )). 



3. CONCLUSIONES 

En el método introducido por Barbosa y Doherty (1987), las variables de 
composición surgen después de analizar la manera en que los potenciales 
químicos se combinan (condición de equilibrio). Por el contrario, en el método de 
los multipicadores de Lagrange y en el de las transformadas de Legendre, la 
elección de las variables que se usarán al describir la composición del sistema se 
hace antes de obtener las expresiones de los potenciales. Esto es lo que permite 
encontrar una expresión sencilla para la función termodinámica que define los 
estados de equilibrio en el caso de reacciones tales como las esterificaciones 
(Pérez Cisneros y otros 1997), es decir cuando los elementos están siempre 
combinados. 

Cuando los "elementos" se eligen de tal manera que la información acerca de la 
composición del sistema no está codificada — es decir, cuando los elementos 
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coinciden con las especies de referencia del método de las combinaciones 
lineales — los resultados convergen y corresponden a la misma transformación de 
Legendre. 

La ventaja de recurrir a las transformaciones de Legendre radica en el hecho de 
que, durante la demostración, el sentido físico de las variables permanece 
explícito. Esto permite hacer una interpretación de la función auxiliar de 
Lagrange, la cual tiene un significado físico claro: es la transformada de 
Legendre completa del potencial de Gibbs respecto de las variables de 
composición de los elementos, y su valor es cero para estados de equilibrio. 

Las variables de composición transformadas son aquellas combinaciones lineales 
independientes de los grados de avance de las distintas reacciones. También son 
las variables conjugadas de los potenciales de los elementos. Dichos potenciales 
son: (1) las variables intensivas que acompañan a ciertas combinaciones lineales 
de los números de moles (apéndice A); (2) los multiplicadores de Lagrange, 
cuando se minimiza el potencial de Gibbs sujeto a restricciones estequiométricas; 
y (3) las variables conjugadas de las variables de composición de los elementos, 
cuando se obtiene la transformada de Legendre completa de G. 
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APÉNDICE A 

VARIABLES DE COMPOSICIÓN TRANSFORMADAS Y GRADOS DE AVANCE 

Tomando como base el sistema estudiado en la sección 2, aquí se demostrará que las 
variables de composición transformadas son las combinaciones lineales de las variables 
de composición independientes de los grados de avance. 

Separando de las cantidades iniciales de B y C las que participan en cada reacción, 
obtenemos 



ft = 



„1 „01 „ „0 „1 „01 

n c -n c n A - n A n B - n B 



-1 
nh = n°c-¿;i; 
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n A = rci + fr; 

ni = ng+tú 



& = 



n r 



„2 _ „02 

n B n B 



<-&, 



n B -n^ 



„2 _ „02 

n c -n c 



n D = n D - £ 



2, 



llr 



n c z + 2{ 2 . 



Y, eligiendo combinaciones lineales adecuadas de las variables de composición, 

n A + n c + 2n B = n A + £. + n°¿ -£i+n c 2 + 2& + 2 (re£ - &); 

= n A + n c + 2n B ; 

n B + n c + n B = n^ 1 + ft + nf - & + n°¿ - ft + n° c 2 + 2 & + ^d - 6; 

= nl + n% + n%. 

A estas variables no corresponde ningún £: no avanzan. 

Adviértase que el punto clave del razonamiento es la separación de las cantidades 
iniciales de B y C en dos partes cuyas extensiones no son arbitrarias. 

APÉNDICE B 

COMPARACIÓN DE DOS TRANSFORMACIONES 

Las transformaciones de las variables de composición en los métodos de los 
multiplicadores de Lagrange y de las transformadas de Legendre es previa a la aplicación 
de esos métodos y puede ser escrita como sigue: 



<n A \ 
n B 
n c 
n-D ) 



n, 



n. 



Por el contrario, la transformación correspondiente a la técnica de las combinaciones 
lineales es el resultado de la aplicación de esa técnica y se puede expresar como sigue: 



n A 

n B 



B 



n c 
n D 



n A 
ñ B 



Los casos en que conducen al mismo resultado pueden ser hallados resolviendo la 
ecuación 
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(n K \ 








n B 


-T A 1- 


(b n b í2 \ 


(n c 


n c 


\n B ) 


U21 b 22 ) 


Ud 


^D , 









Las componentes de C son: 



Cu Ci 2 C13 C14 
C21 C 22 C23 C 2 4 



1 -6n -6 12 
1 -621 -622 



Por lo tanto, las transformaciones llevan al mismo resultado cuando la matriz C puede 
ser dividida en una matriz diagonal y otra que es la opuesta de B, es decir, cuando los 
elementos se toman del conjunto de los componentes. 

APÉNDICE C 

FORMAS DIFERENCIAL E INTEGRAL EN LA TÉCNICA DE LAS 
TRANSFORMACIONES DE LEGENDRE 

En la sección 2, el problema fue resuelto por medio del uso de transformaciones de 
Legendre, usando ecuaciones diferenciales. El mismo problema puede ser resuelto 
usando esas ecuaciones en su forma integral. En este caso, el método viene a ser el 
recíproco del introducido por Barbosa y Doherty (1987): dadas las ecuaciones para 
obtener n* a partir de n, el objetivo es hallar la relación entre C, y ¡j, . 



G 



(& 6) 



10 12 
111 



n B 

V n D ) 



= 



(Ma~U /¿b-£b ^c-^a-^b ¿Ud-2£ a -£b) 



n B 

n c 

V tid ) 



Esta condición debe ser cumplida cuando no todos los números de moles son nulos 
(solución trivial), 

( ¡x K fi B nc Ht> ) = ( U £b U + £b 2 £ A + ¿f B )• 
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